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Einstein-Smoluchowski diffusion, damped harmonic oscillations, and spatial decoherence are spe¬ 
cial cases of an elegant class of Markovian quantum Brownian motion models that is invariant under 
linear symplectic transformations. Here we prove that for each member of this class there is a pre¬ 
ferred timescale such that the dynamics, considered stroboscopically, can be rewritten exactly as 
unitary evolution interrupted periodically by an entanglement-breaking measurement with respect 
to a fixed overcomplete set of pure Gaussian states. This is relevant to the continuing search for the 
best way to describe pointer states and pure decoherence in systems with continuous variables, and 
gives a concrete sense in which the decoherence can be said to arise from a complete measurement 
of the system by its environment. We also extend some of the results of Diosi and Kiefer to the 
symplectic covariant formalism and compare them with the preferred timescales and Gaussian states 
associated with the POVM form. 


Although it has been widely studied for the better 
part of a century, the dynamical equations for Marko¬ 
vian quantum Brownian motion (QBM) [IH3] were not 
solved in full symplectic generality until relatively re¬ 
cently [IHB]. By symplectic generality, we mean a class of 
dynamical equations for a quantum state that is invariant 
under linear symplectic transformations of phase space, 
i.e., transformations of {x,p) that are linear and pre¬ 
serve the symplectic form, in close analogy with Lorentz 
covariance. Equation 0 below describes the minimal 
such class of dynamics that subsumes the harmonic os¬ 
cillator, frictionless spatial decoherence, and Einstein- 
Smoluchowski (frictionful, noninertial) diffusion. QBM 
models have been an essential testbed for understand¬ 
ing decoherence and the quantum-classical transition in 
systems with continuous degrees of freedom [THU] , es¬ 
pecially in the special case that can be generated by an 
environmental bath of oscillators coupled linearly to po¬ 
sition. 

In discrete systems, pure decoherence [m US] serves 
as a platonic ideal to which many real-world systems are 
close approximations. Pure decoherence of a system oc¬ 
curs when the system’s reduced dynamics take the form 
of a dephasing channel with respect to some orthonormal 
pointer basis HU El [HI- However, the natural analogs 
of pure decoherence and the pointer basis remain elusive 
for systems with continuous degrees of freedom; several 
different definitions for the pointer basis in such cases 
have been proposed US El [MM] but none are widely 
accepted. 

It is now well appreciated that complete pure deco¬ 
herence in discrete systems is equivalent to a complete 
measurement of the discrete variable by the environment. 
In this article we generalize this to continuous systems. 
We prove that time-homogeneous Markovian QBM, when 
considered as a stroboscopic evolution between discrete 
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times, is exactly equivalent to unitary evolution punc¬ 
tuated periodically by a fixed positive operator-valued 
measurement (POVM) with respect to an overcomplete 
set of Gaussian wavepacket states. In other words, the 
non-unitary component of QBM evolution is described 
by an entanglement-breaking [251126j Gaussian measure¬ 
ment carried out by the environment. This result can be 
straightforwardly extended to the case where the Hamil¬ 
tonian and Lindblad operators are time-dependent [27] . 
This immediately recovers earlier work showing that, un¬ 
der QBM, the Wigner function becomes strictly and per¬ 
manently positive in finite time for an arbitrary initial 
state |4l|2Hl|29]. 

It is tempting to suggest that the Gaussian wavepack- 
ets composing that POVM should reasonably be iden¬ 
tified as the precise pointer states of this open system 
evolution. However, even when the dynamics are cast 
into the POVM form, there is still significant freedom 
to choose the preferred pointer states. The choices are 
close to, but generally distinct from, the (also Gaussian) 
pointer states suggested by Diosi and Kiefer dUHl] and 
by Zurek and collaborators m E] . Such alternatives 
are related to different criteria of classicality that have 
been studied in the past, such as the Wigner function 
becoming positive [SI1H1121], the Glauber P function be¬ 
coming positive [281130] (and no more singular than a <5 
function |3I]), or that the quantum state p is expressible 
as an incoherent mixture of Gaussian states [31] . Thus, 
our results add to the zoo of possible classicality crite¬ 
ria and associated pointer states used to understand de¬ 
coherence in QBM, although most sensible choices are 
closely related for dimensional reasons. 

In Section I we state the main result after introducing 
the minimum necessary notation, and in Section H we 
give a proof. In Section HI we connect this to other 
preferred states and timescales discussed in the literature, 
especially to the work of Diosi and Kiefer. In Section IV 
we offer concluding discussion. A summary of symplectic 
QBM using our notation can be found in Appendix A, 
including a description of important special cases. In 
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Appendix B we explicitly compute the POVM form of 
the dynamics for the common special case of a lightly 
damped harmonic oscillator. 


I. MAIN RESULT 

Markovian QBM in the symplectic general form is 
given by the master equation dtpt = >C[pt] where the man¬ 
ifestly covariant Linblad superoperator is 

^ [d“,[d^pt]], ( 1 ) 

for the density matrix pt of a single continuous degree 
of freedom m- Here, a = {x,p) is a vector operator 
for a point in phase space, repeated indices are summed 
over the two phase-space directions (a = x,p), Dab is 
a positive semidefinite 2 -by -2 matrix with real entries, 
and Kab is a 2 -by -2 matrix with real entries satisfying 
(K^)^ < 2DabD°‘^. Indices are raised and lowered with 
the symplectic form using the anti-symmetric Levi-Civita 
tensor e“^, e.g., = e°‘‘^Kcb- (They behave just like 

Weyl spinors.) 

We set h = 1 and introduce notation where 
and are replaced by boldface K, D, and a. We de¬ 
fine 7 E -TrK/2 = -K\I2 = A:„“/2 and H e K + 7 I 
so that H is traceless and 0 < 7 ^ < |D|, where I-] de¬ 
notes the determinant. Let V be an arbitrary matrix 
with unit determinant representing a canonical linear 
transformation, i.e., a (classical) linear transformations 
on phase space which preserves the symplectic form. We 
denote the associated quantum unitary evolution by LV, 
so that {ul^aU'v) = V {a) and = Uv-i. The corre¬ 
sponding superoperator on the space of density matrices 
is Hv[p] = U\jpij\^. Let |a) represent the normalized co¬ 
herent states with phase-space mean a. = {x,p), so that 
|V;a} E Uv\V-^a) are the pure Gaussian states [35] 
parametrized by their mean a = (a) and their 2 -by -2 
covariance matrix V where V e VV^/2. 

Lastly, for s > 0, 

/^(V;a|p|V;a>.|V;sa>(V;sa| (2) 

represents the entanglement-breaking channel Q given by 
the Krauss operators |V;sa)(V;a|. This is a 

POVM measurement with respect to the overcomplete 
basis |V;q:] - formally, a frame |3dj - where the mea¬ 
surement outcome o: is followed by a preparation of the 
state |V;sa](V; sq:| (and then forgotten). A minimally 


^ A channel on a Hilbert space A is defined to be entan¬ 
glement breaking if it always produces separable states when 
operating on one part of an entangled pair, that is, when 
"S*is separable for all '^aB' This is true if 
and only if all its Krauss operators are unit rank m- 


disturbing measurement corresponds to s = 1, while dila¬ 
tions and contractions of phase space are given by s > 1 
and s < 1, respectively. 

Theorem. There is a characteristic time T > 0 defined 
as the unique positive solution to ICrl = (1 + e"^'’'^)^/4, 
where 

CtE r‘e^*^De^^"dr. (3) 

Jo 

Moreover, integrating the dynamics ([T]) forward by T in¬ 
duces a completely positive (CP) trace-preserving map 
that sequentially evolves, measures, and prepares the sys¬ 
tem, 

= (4) 

where o denotes composition, and where V = Ct/( 1 + 
e“^^^), s = and R = 

Remark. This can be rewritten in several ways: 

= Hv o (hj o Hv-ir[p] 

= nn o (5) 

= f^(RFW-,(x\p\RFW-a) 

X IV; sRa](V; sRq:| , 

or generally as 

[p] = Hri-t- o o Hrt-[ p]. (6) 

where 

V - r-A-ATH Ct 

= r-(i-Uv(RT)-(i-’') (7) 

= R-(i-^)V 


for any r € [0,1], with special cases Vg = R’^V and 
Vi = V. Since the dynamics are Markovian, this can be 
iterated, e.g., 

e2T£ = $AoHRO$^oHR, (8) 

and so on for with n any positive integer. In other 
words, quantum Brownian motion can be understood 
stroboscopically as an iterated phase-space measurement 
interspersed with unitary evolution. 


II. PROOF 


Under Q, the Wigner function Wp corresponding to 
the state pt is known to obey a Klein-Kramers dynam¬ 
ical equation EHSj, which can be written in sympectic 
covariant form as 


dtWpia) = 


-K\daa‘’ + -D^^dadb 


Wp{a), (9) 
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where da = dlda°‘. Below we will work within the space 
of functions over phase space using the convolution oper¬ 
ator * and the composition operator o. In this context, 
matrices are taken to represent the functions obtained by 
matrix multiplication with the phase space point a, i.e., 
V(a) = Va and {Wp o V)(a) = Wp{'VoL)_ 


In this notation, the exact solution to (24) for any ini¬ 


tial Wigner distribution Wp{a) is known to be nn 
where C( is given by @ and 


( 10 ) 


GB(a)E 




27rv^ 


( 11 ) 


is a normalized Gaussian smoothing kernel for any pos¬ 
itive semidefinite covariance matrix B. Note the impor¬ 
tant special case of unitary evolution, 7 = 0 = D, for 
which kkgt£[p] = Wp o 

The key idea in the proof is that the POVM-and- 
prepare channel $1 with respect to the coherent states 

|a) = |I;q:) corresponds, in the Wigner representation, 
to a convolution with the kernel = Gj. That is, 
ll^$l[p] = G 2 J * Wp. More generally, dissipation in phase 
space can be accounted for by considering the modified 
channel 

= /^ka)(a|p|a)(sa| ( 12 ) 


and calculating the corresponding Husimi Q function 

Q$£[p](a) = ^ (a|^f[p]|a) 

rip , (13) 

= J £QpmHsf3)f 

where we have changed integration variables to /x = 
s/3. Using the fact that |(q:|/ 2 )|^ = exp[-(Q; -/x)^/2] = 
27rGi(Q: - fi) we obtain 

Q$»[p] = G 21 * (s"^Qp ° s~^I) (14) 


The Husimi Q function is just the Wigner function 
smoothed by a Gaussian, Q = Gj* W. Using these di¬ 


rectly checkable identities 

(Xx-y)oA=|A|-(XoA)x-(yoA), (15) 

Gb ° A ^ = |A|-Gabat, (16) 

Ga*Gb = Ga+b, (17) 

we can deconvolve both sides of ( [l4| ) to get 

ll^^£[p]=G(i,,.)j.(s-^WpOs-il) (18) 


This can be generalized to a POVM of Gaussian states 
|V;q;) = {7v |l; with any linear symplectic trans¬ 

formation V by first applying an appropriate unitary 
Hv-i, executing the measurement, and then applying the 
inverse unitary Hv: 4>^[p] = Hv o o Hv-i [p]- In the 
Wigner representation this is 


ll^^iL[p] = ll^nvo$£onv-i[p] 


4))oV- 


(19) 


= [p] ° V 

= (G'(i+s2)l"-(s"^14^n^-i[p]°s 

= (%.«=)! °V-i).(s-X°«-'l) 

Augmenting this with unitary evolution Ur, corre¬ 
sponding to an arbitrary linear symplectic transforma¬ 
tion R gives 

Il^$AonR,[p] = G^^^g 2 )v * (s Wp o s R ) (20) 


We can then reproduce the solution (10 1 by choosing 
V = Ct/(1 + s^), s = and R = . However, 

this is only possible when ICtI = (1 + s^)^/4, since we 
have assumed V is a linear symplectic transformation (so 
|v| = |V| /4 = 1/4). Let us prove that this requirement 
uniquely determines T so long as 7 is not extremal (i.e., 
so long as 7 ^ < |D| rather than 7 ^ = |D|). 

First, note that since D is positive semidefinite, the 
integrand of (|^ is also positive semidefinite by con¬ 
struction. By the Minkowski determinant theorem (see, 
for example. Ref. [M]), the determinant function obeys 
\/iArU^ > -y/jAl + \/j^ for positive semidefinite 2 -by -2 
matrices, so 


|Ct|> 


f y/|e*^De*^ 

Jo 


|D| 


( 21 ) 


= (l_e-27*)2^^ 


M 

472 ’ 


and likewise 


at 


( 22 ) 


From this one can show that 4|Ci|/(l + e 27*^2 ijj^eases 
monotonically with t and moreover 


4|Ct 


(1 + e-27‘)2 


> tanh^(7t) LJ. 


(23) 


For 7 0, the function tanh^( 7 t) starts at zero and ap¬ 

proaches unity as t ->• 00 . Thus for |D| 7 0, we have that 
4 |Ct|/(1 + 6 “^'’'^)^ = 1 for some unique finite T, except 
for the extremal cases 7 = ±-y|D| (including D = 0) for 
which r -> 00 . 
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III. OTHER STATES AND TIMESCALES 


In this section, we consider the preferred timescales 
and pointer states discussed by Diosi and Kiefer 
IdOj and others PEHlCini in the context of a frequently 
studied special case of QBM, and generalize them to the 
symplectic covariant formalism. We compare them to 
the preferred states and timescales associated with the 
POVM form for the dynamics derived above. We do not 
necessarily expect closed-form expressions for arbitrary 
K and D, but we can nevertheless show that the preferred 
quantities are well-dehned and generally distinct. 

The frequently studied special case of the dynamics 
can be described as momentum diffusion and spatial de¬ 
coherence that is frictionless (7 = 0 ) and spatially homo¬ 
geneous. This is often called simply quantum Brownian 
motion, but we will call it pure spatial decoherence to 
distinguish it from the general case, Q. It is defined by 
setting = I/m, = D, and all other coefficients of 
and D°'^ to zero. The dynamical equation for the 
Wigner function reduces to 


dtWt{x,p) = 


dx + „ dp 
m 2 . 


Wt{x,p). (24) 


We can obtain the dynamical equations for Qj by 
making the replacement Wp Qj and D ^ in 
where 

d|,e D-(KY +YK'^). (29) 


Likewise is true for , by making the replacement Wp 
and D ^ in Q, where 



D^eD + (KY + YKT). 

(30) 

The solutions are 




(31) 



) (32) 

where 

EFEe‘^Ye‘^^-Y 

= [KY + YK'^] dr 

(33) 

These 

can be checked using (27-28|), (|15 -I7), 

and ( 10 ). 


Pure spatial decoherence is often obtained mathemat¬ 
ically from an explicit model of the environment as a 
thermal bath of oscillators coupled linearly in x, followed 
by taking the large-temperature limit ISHUE- It also 
well describes the dynamics taken by a test mass sub¬ 
jected to collisional decoherence [SI EH |361 El] from an 
environment of lighter particles |38j , blackbody radiation 
[5], or low-mass dark matter [55] , 

In order to extend the results of Diosi and Kiefer to 
symplectic generality, we recall that the Husimi Q func¬ 
tion and (when it is well-defined) the Glauber P function 
can, for any state |'(/>), be usefully generalized [40] to 

^ (25) 

p= j da P''^\a.)\'il},OL){-il},a\, (26) 

where |^/',Q:) = Ta is the state I'i/) translated in phase 
space by a. These reduce to the original Q and P func¬ 
tions when |' 0 ) is the coherent state centered at the ori¬ 
gin in phase space, i.e., the ground state of the harmonic 
oscillator, |I;q: = O). In the context of quadratic Hamil¬ 
tonians, it is natural to concentrate on the case of gen¬ 
eralized Q and P functions for which l-ijj) = |Y;a = O), 
i.e., a general Gaussian state centered at the origin with 
covariance matrix Y, and adopt the shorthand notation 

QJ = and P^ = These can be related 

(see Appendix |A 4| ) to the Wigner function by 

Qj=G^*Wp, (27) 

Wp = G^*Pj (28) 


A. Unraveling the Glauber P function 


As observed in Refs. EEO], is noteworthy be¬ 
cause, when Pjtc[p^ {oP) is defined, the system is described 
by a (classical) probability distribution diffusing over a 
set of pure Gaussian states with a preferred covariance 
matrix Y. However, this interpretation is only viable 
when the diffusion matrix e D + (KY + YK^) for 

the Glauber function P^ is positive semidefinite. This 
restriction > 0 dehnes a region in the space of possi¬ 
ble Gaussian pure-state covariance matrices Y. 

Note that this region may be empty for some choices 
of dynamical parameters]^ For instance, if = P = 
for real p and all other parameters vanish, then 
no choice of Y makes positive. In this case, the 
dynamics continuously squeeze phase space toward one 
axis, so that an initial Gaussian state becomes arbitrarily 
squeezed with increasing time, and hence eventually not 
expressible as a mixture of Gaussians with fixed, finite 
covariance matrix Y. 

If there is a region of compatible Y with finite vol¬ 
ume, it is then possible to look for an additional crite¬ 
rion that would prefer some states in this region over 


^ We observed numerical and analytic evidence that the region 
DA > 0 has strictly positive volume for all but a measure zero 
subset of the dynamical parameter space. That is, it appears 
that there is always a choice of Y that allows for the diffusive 
interpretation with DT for almost any dynamical parameters. 
However we could not find a proof of this. 
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others. One choice is to find the pure initial state \4’) 
that minimizes the instantaneous linear entropy produc¬ 
tion dS'L/dt, where = 1 -Trp^. This is motivated 
by the intuitive notion of the predictability sieve; pointer 
states are the quantum states that are most stable un¬ 
der interactions with the environment [H [m Eg, and 
thereby produce little entanglement entropy. 

In general, the linear entropy production is minimized 
by choosing the most squeezed dimension of the initial 
state to be along the largest eigenvector of the diffusion 
matrix D. We follow Kiefer et al m and Diosi and Kiefer 
m and parametrize the most general unit-determinant, 
positive semidefinite matrix Y according to the complex 
number = 17^ + i^i (with > 0): 

This corresponds to a spatial wavefunction (a;|Y ;a = O) = 
(17^/271)^'^^ exp(-f7a;^/4). In the special case of pure spa¬ 
tial decoherence, the linear entropy production is pro- 


I7j 


V Dm 



17r 
V Dm 


FIG. 1. (Color online) Some candidate Gaussian 
pointer states with spatial wavefunction {x\Y ; q : = 0^ oc 
exp[-(S7R -I-if2i)a:^/4] and corresponding covariance matrix 
given by ( |34[ ). Wavefunction normalizability requires that 
S7 r > 0. The arc corresponds to the states Y = Vr compat¬ 
ible with a POVM form, Q, for the dynamics, (24|. It is 
parametrized by r s [0,1], running down from (a) Vi = V to 
(b) Vi /2 to (c) Vo = R ^V. Under infinitesimal evolution, 
the state that can be best approximated by a pure-state ac¬ 
cording to the Hilbert-Schmidt norm is labeled by (d). States 
within the shaded region satisfy > 0, allowing for an in¬ 
terpretation in terms of a diffusing probability distribution 
over pointer states. The state producing the minimal lin¬ 
ear entropy (equivalently, maximum IIr) from that region is 
labeled by (e). 


portional to I7 r, and one can check for pure spatial 
decoherence ( p4| ) that I7 r is maximized by the choice 
17 = - i)\/Dm under the constraint that is 

positive El- 

One may alternatively consider the pointer state se¬ 
lected by the principle of Hilbert-Schmidt robustness 
[niiiiiEniiinii]- The time-dependent pure states 
that best approximate, according to the Hilbert-Schmidt 
norm, the impure state pt evolving under QBM can be 
shown to solve [30] 


^\^t) = [c[\^t){H] - (V't|-C[|V'i)(V'd]lV't))lV't)- (35) 


For pure spatial decoherence (241, the unique stationary 
solution to this non-linear equation are candidate pointer 
states [TTl 1331130113311331130] . They are all equivalent up 
to translations in phase space, being given by a Gaussian 
wavepacket with a covariance matrix specified by 17 = 


{l-i)V^ m El- 

In Fig. these two preferred states and the preferred 
region associated with the condition > 0 are com¬ 
pared with the one-parameter family of Gaussian states 
characterized by the covariance matrix given by Q 
that are compatible with the general POVM form ([^ for 
pure spatial decoherence. 


B. Positivity times 


Diosi and Kiefer also calculated [28] the characteristic 
times T\y and Tp at which the Wigner function Wp and 

the traditional Glauber function Pp = of an arbitrary 
quantum state became strictly positive under pure spa¬ 
tial decoherence (24). As they conjectured was possible, 
this was extended to symplectically general QBM dy¬ 
namics by Brodier and Ozorio de Aleida [3| . The Wigner 
positivity time Tyy is strictly a property of the dynamics 
in the sense that the time is independent of the initial 
Wigner function (so long as it is pure and not already 
positive) [3]. Once a Wigner function is positive, it re¬ 
mains so indefinitely under QBM dynamics. Here we col¬ 
lect these results and likewise treat the positivity of the 
Glauber Pj function, corresponding to Gaussian kernels 


with arbitrary covariance matrix Y, in symplectic gener¬ 
ality. 

First note that 


W,cT[p]=Gc,*{e^^'Wpoe-^^) 

= {{Gc,°e^^)*Wp)oe-^^ (36) 


The key idea is that convolving the Wigner function by 
a Gaussian Gj yields the Husimi Q function, and the lat¬ 
ter is always positive |28| , but convolving with a sharper 
Gaussian (e.g., G 1 / 4 ) will never produce a positive func¬ 
tion from a nonpositive pure state (or vice versa) [3|. The 
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classical flow and the multiplicative factor > 0 
do not change the positivity of a Wigner function, so con¬ 
volution by Gct will make a Wigner function for a pure 
state positive if and only if |Ct| > |l| = 1/4. 

Thus by examining (36), Tw is defined by 


|Ct„ 


\=-e 

4 


-i'yTw 


(37) 


Using (21), once can show that is finite (except for 


the extremal case 7 = -y/|D|). Likewise, the generalized 
Glauber P function associated with a set of preferred 
Gaussian state with covariance matrix Y is given (when 
it exists) by Wp = * Pj 


p ’ 


so 




= I 




(38) 


(^e-*K(Ct-Y)e-*KT * fUp j o e 

It is guaranteed to exist and be positive at the charac¬ 
teristic time Tp defined by 


I 14 


(39) 


Different choices of the preferred covariance matrix Y 
lead to different times upon which the associated Glauber 
function becomes positive. For pure spatial decoher¬ 
ence (24), Tw = = T/\/2 [2S], and Tp is the 

solution to 12 m^ + GmUit + (Oj + 0 ^)t^ - = 0. 

Many other extensions are possible. For instance, it 
is clear that one could calculate the time at which the 
Gahill R function, which continuously interpolates be¬ 
tween the P, W, and Q function [47], becomes positive 
for different values of the interpolation parameter. Like¬ 
wise one could define pointer states of the system to be 
the Gaussians with covariance rnatrix Y such that diffu¬ 
sion in the Glauber function P^ is preferred according 
to a criterion other than minimizing the linear entropy 
production. None of these stand out as definitive notions 
of classicality. 

It is worth emphasizing that the condition Wp > 0 is a 
property strictly of the quantum state p itself, whereas 
the condition Pj > 0 (and most other criteria based 
on the functions P^ or Q^) are dependent on the di¬ 
mensionful choice Y. The traditional definitions for the 
Glauber and Husimi functions Pp = Pj and Qp = Qp do 
not avoid this because they depend implicitly on a length 
scale Gx used to define the identity matrix, 


1/ 

' 1 



0 \_j 

I 

2 ' 

.0 


[ 0 


[ 0 


0 


(40) 


where = {x^) - (x)^ = {a = 0|a:^|Q: = 0) is the spa¬ 
tial variance of the coherent states. This length scale is 
usually taken from the dynamics (most often, the width 
of the oscillator potential) and is not a property of the 
state p alone. Likewise, it is often taken for granted that 
squeezed states are nonclassical (e.g., |15|), but squeezing 
is always relative to an assumed scale separate from p. 


IV. DISCUSSION 

General QBM dynamics arise from the lowest order 
terms in the Taylor series expansion of a smooth Hamil¬ 
tonian for a Markovian open system, giving them similar 
conceptual importance and pedagogical usefulness as the 
harmonic oscillator has in the study of closed quantum 
system. The linearity of the Lindblad operators means 
their influence can be described as continuous weak moni¬ 
toring of the phase-space variable a = {x,p) OHHHH]. In 
this work we have shown how these weak measurements 
add up to a single strong measurement, on a timescale 
that characterizes the dynamics, in the form of a POVM- 
and-prepare (entanglement-breaking) quantum channel. 
The symplectic generality exhibited here will be very im¬ 
portant for extending, to Markovian open systems, exist¬ 
ing quasiclassicality theorems [311 1501153] that apply to 
evolution generated by any closed-system Hamiltonian 
that is sufficiently smooth to be treated as approximately 
locally quadratic, rather than just for representative toy 
environments like baths of harmonic oscillators. 

The various forms exhibited in ([^ suggest that the 
overcomplete set of Gaussian states forming the POVM 
are only defined up to a sort of gauge freedom 

(0,{|V;a)})-(t,{|R=^V;a)}) (41) 

for any t € [0,T]. On the other hand, the time scale 
T associated with the dynamics is independent of this 
freedom 0 

This is the continuum analog to a preferred basis am¬ 
biguity that can be found in the simpler case where the 
dynamics are described stroboscopically by unitary evo¬ 
lution punctuated by a simple projective measurement 
of a discrete variable. See Fig. [^ The ambiguity arises 
because of the periodic measurement events and the non¬ 
trivial unitary evolution of the system in between them; 
it is not a special property of continuous system. 

Given this freedom, as well as the alternative pointer 
state criteria discussed in the previous section, it is not 
clear whether the pointer states of quantum Brownian 
motion are best understood in terms of the POVM form 
for the dynamics presented above. However, this basis 
ambiguity may play a conceptual role in any future sat¬ 
isfactory notion of pointer states when decoherence is 
taking place alongside unitary evolution. In other words, 
one should be suspicious of the intuition that there is 
a single true preferred basis (whether overcomplete or 
otherwise) that one might develop from studying simple 
models of pure decoherence in an orthonormal basis. 


® Of course, all POVMs with respect to the frame of Gaussian 
states with a covariance matrix Y are equivalent if we allow 
them to be supplemented with an arbitrary unitary immediately 
before and after the measurement. The one-parameter family in 
is notable because it is merely inserted at some point in the 
normal unitary component of the evolution. 
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FIG. 2. (Color online) The real part of the Hilbert space of a two-dimensional system. The system evolves stroboscopically 
according to a unitary Uo which is periodically interrupted every time step T by a projective measurement in the basis 
{|Si), IS' 2 )}. After two measurements, there are four possible outcomes: (1,1), (1,2), (2,1), and (2,2), corresponding to the 
final states \vij) = |5jX<S'j| Uo |5'iX'S'i| Uo |no). The basis {|5i)} plays the role of a pointer basis, but this evolution is equivalent 
to periodic measurements in the basis I'S'i)} with a constant time offset t e [0,T] (whether or not the outcomes i are 

remembered or forgotten). In this sense, the pointer basis for this system is only sensibly specified relative to an offset time, 
i.e., as a pair \Si)}). 


Although the characteristic time T of the POVM de¬ 
scription differs (by a factor of order unity) from the ex¬ 
act times Tw and Tp at which the Wigner and Glauber 
functions become positive, the positivity manifestly im¬ 
plied by the former description is arguably more trans¬ 
parent; the Wigner function following a POVM-and- 
prepare channel in a Gaussian state basis is obviously 
positive. It also suggests new approaches to understand¬ 
ing phase-space positivity, or the quantum-classical de¬ 
scription more generally, in continuous-variable systems 
that aren’t described by QBM. 

The symplectic covariance of our results offers illu¬ 
minating generality compared to earlier discussion of 
pointer states in special cases of QBM [HHUl [H HI] ■ For 
example, the pointer states associated with the POVM 
form for the damped harmonic oscillator dynamics are 
a generalization, to arbitrary damping, of the coherent 
states that were identified as pointer states in the under¬ 
damped limit (7 « uj) [T4] . 

It is notable that, except in the extremal case 7 ^ = |D|, 
the preferred states associated with the POVM form 
for the QBM dynamics are always well defined and un¬ 
ambiguous, a result which may also apply to pointer 
states associated with Hilbert-Schmidt robustness. In 
contrast, the predictability sieve often produces singu¬ 
lar pointer states like the position eigenstates with di¬ 
vergent momentum dispersion HD: unless supplemented 
with additional cumbersome principles such as a Glauber 
P function dispersion interpretation, or a finite-time av¬ 
eraging scheme. Of course, no elegant principle exists 
that unambiguously identifies sensible pointer states (or 
their nonexistence) for arbitrary dynamics, and the pre¬ 
dictability sieve appears to offer more guidance there. 


Our most restrictive assumption has been that the dy¬ 
namics are Markovian and time homogeneous. One way 
to relax this is by allowing D and K to vary with time, 
possibly in a way that depends on the initial state. This 
will lead to straightforward modifications of C*, T, and 
V, and it is still possible to describe the dynamics strobo¬ 
scopically as a Gaussian-state POVM-and-prepare chan¬ 
nel [57]. In this case the preferred states and timescales 
are generally not determined solely by the dynamics, but 
also by the initial state. 

This might be extended to cover more general models 
of non-Markovian dynamics, like the finite-temperature 
bath of linearly coupled oscillators of Galdeira-Leggett 
[aia 112] . However, in such cases D is not necessarily 
positive-definite, and this can interfere with construct¬ 
ing the POVM form. Indeed, this form must breakdown 
when the memory of the environment exceeds the strobo¬ 
scopic time interval. (Trivially, a finite bath has a finite 
global recurrence time and so will eventually restore any 
initial non-classical superposition states of the system.) 

It would be especially interesting to see if pointer states 
can be identified in non-Markovian dynamics such that 
the environment’s memory consists only of the classi¬ 
cal history of those preferred state, e.g., if the evolution 
can be described as an iterated sequence of POVM-and- 
prepare channels depending on previous outcomes. On 
the other hand, it is hard to see how pointer states could 
be usefully defined for the more general non-Markovian 
case where the coherence information between pointer 
states feeds back from the environment into the system; 
in that case, one would rather say that decoherence had 
not been effective and there simply are no pointer states. 
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Appendix A: Symplectic QBM 


In this appendix we briefly review QBM along the lines 
of Brodier and Ozorio de Almeida [4] and Robert [6] , but 
in context and notation better suited for our purposes 
that especially emphasizes manifest symplectic covari¬ 
ance. See also the appendix to Ref. [5S], Chapter 6 of 
Ref. [5], and a forthcoming pedagogical treatment [27] . 


1. Linear symplectic transformations 


SL 2 (K) down to a two-dimensional manifold that can be 
parametrized as in (34). 

(Note that in the context of Gaussian quantum phase- 
space distributions, some authors differ from our conven¬ 
tion by setting h = 2. In this case, the coherent state 
satisfies {x^){jP) = cr^cTp = ft^/4 = 1, i.e., the covariance 
matrix is I rather than I = 1/2 [55]. 1 


2. Quantum Brownian motion 

A single quantum continuous degree of freedom un¬ 
dergoing open-system, time-homogeneous, and Marko¬ 
vian dynamics, forms a quantum dynamical semigroup 
[SUIT] described by a Lindblad master equation 


For a single continuous classical degree of freedom, 
phase space is a two-dimensional vector space equipped 
with the symplectic form as represented by the antisym¬ 
metric Levi-Civita symbol e“^ (with = epx = +1 = 
= -exp)- Analogously to Lorentz indices, sym¬ 
plectic indices are raised and lowered by contracting 
with the second index in the symplectic form. Unlike 
for a Lorentzian metric, which is symmetric, the anti¬ 
symmetry of the Levi-Civita symbol means there is an 
overall sign flip depending on which index is raised and 
which is lowered in a contraction: XaY"' = eabX^Y°‘ = 
-ebaX^Y°' = -X^Yb- The symplectic indices behave 
just like Weyl spinor indices, which are reviewed in 
many introductory quantum field theory textbooks (e.g.. 
Ref [M]~). 

Classical Hamiltonian evolution is given by a time- 
parametrized family of symplectomorphisms (canonical 
transformations) on phase space, which are characterized 
by the fact that the Jacobian of the transformation pre¬ 
serves the symplectic form at each point. We concentrate 
on the local dynamics of smooth Hamiltonians, so we are 
most interested in the symplectic linear transformations, 
which, for one degree of freedom, are effected with the 
Lie group of 2-by-2 real matrices 1"“^ with unit deter¬ 
minant, SL 2 (K). These preserve the symplectic form: 
^ab ^ equivalently, e = YeYU 

Each (classical) symplectic linear transformation Y is 
associated with a corresponding quantum unitary trans¬ 
formation U-Y = exp(^-iZab&°'Oi^l‘2), where Z is the ma¬ 
trix Z\ = e°‘'^Zcb and Y = e^. In particular, [U^aUY) = 

Y {a) and LV-i = Uy = When this unitary acts on 
a coherent states |a) = |I/2 ;q;), it transforms it [5^ [55] 
into a (generally squeezed) Gaussian state centered on 
the classically shifted point in phase space: Uy |q:) = 
|Y;Yq;), where Y e YY^/2. Here, the pure Gaussian 
states |Y;q:) are parametrized by their mean a = (a) 
and their 2 -by -2 positive definite, unit-determinant co- 

variance matrix Y with elements Y = /2. 

(Y determines |Y;q;) up to a phase, which is suffi¬ 
cient for our purposes.) The many-to-one mapping Y -> 

Y = YY^/2 collapses the three dimensional Lie group 


dtpt = -i[H,pt] + Y, 




(Al) 


Ideal QBM is the special case when the Hamiltonian is 
quadratic and the Lindblad operators are linear with 
the phase-space operators: H = )^HabcUcU and = 

La^a°', where Hab is a real symmetric matrix and the 
La ^ are complex vectors. (Linear terms in the Hamil¬ 
tonian either can be handed explicitly separately [5] 
or, so long as the quadratic terms Hab are nonzero, 
can be eliminate through a phase-space translation 
{x,p) (x + a;o,p + po)-) We can then change variables 
to Dab = ReZ,(LW)*L(*\ 7 = )* 4 ‘' = 

and Kab = Hab + eabl- (Note that 
our convention for the matrix Dab agrees with Diosi and 
Kiefer mm, but differs by a factor of two from Isar et 
ah, [3|, Dekker and Valsakumar [5S], and others.) 

In the traditional Hilbert space representation this 
yields 

= -^^Kab [d“, Pt]] - \Dab [d“, [d^ Pt]] , (A2) 
Using the Wigner function 

Wp{x,p) = ^ f dAxe~’‘^^^ {x + Ax/2\p\x- Ax/2) 

2t: J 

(A3) 


we can instead express these dynamics in the Wigner 
representation as a Fokker-Planck equation 


dtWp{a) = 


-K\daa^ + -D<^'‘dadb 


Wp{a). (A4) 


Because the equation for the classical phase-space prob¬ 
ability distribution under ideal Brownian motion is iden¬ 
tical to (A4), one can directly read off an interpretation 
of the coefficients. In the absence of diffusion (D = 0), 
the classical equations of motion are 6l = Kct, with 
the classical (possibly dissipative) flow and 
the Hamiltonian component. The rate of dissipation 
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is 7 , where 7 < 0 implies that the environment is pump¬ 
ing energy into the system. The (strictly quantum) con¬ 
straint |D| > ?i^ 7 ^ ensures that diffusion is always suffi¬ 
ciently strong to prevent the dissapative contraction e~*'^ 
from producing violations of the uncertainty principle. 

We briefly mention the most important special cases. 
(For others, see especially Ref. [3].) Normal friction 
(drag) on a free particle of the form x = -{Xlm)x is ob¬ 
tained with RPp = -Xlm < 0 , K^p = 1 /m > 0 , and all 
other parameters zero. These dynamics are dissipative 
since 7 = -K'^^ = X/m > 0. A harmonic oscillator po¬ 
tential is represented by = -muP, and the degree of 
damping (under- or over-damped) is controlled by A. 


Einstein-Smoluchowski diffusion is obtained when nor¬ 
mal friction is supplemented with momentum diffusion 
(£)PP = 2Dq), since the motion is caused by many small 
momentum transfers from molecular collisions. However, 
it takes place in the noninterial limit where the relaxation 
timescale 7 “^ = m/X is short compared to the timescale 
on which observations are made. The fast relaxation 
means the momentum cannot grow, so on large scales 
the position acts as a random walk, oc T. Contrast 
this with the frictionless momentum diffusion dynamics 
(pure spatial decoherence) described by (24), for which 
the position variance (x^) grows like rather than T. 


4. Generalized phase-space distributions 

In the case of a preferred covariance matrix Y, the 
generalized Husimi Q and Glauber P functions are 

Qj{a)^^{Y;a\p\Y;a), (A7) 

p= Jda pJ{a)\Y-cx){Y-a\ (A8) 

These can be related to the Wigner function by first con¬ 
sidering an arbitrary pure state p = |x)(x| and calculating 

QM(x|(“) = ^l(Y;ak>r 

= ^|(Y-ia|l7;|x)r 
= J d/3 W|Y-iaXY-ia|(/3)W^aY-i[|x)(x|](/^) 

= j^d/3Gr(Y-ia-/3)H|^)(^|(Y/3) 

= f d/xGj(Y-i(a-/x))H|^)(^|(/x) 

(A9) 

where we have made use of the fact that Tr[/ 9 p'] = 
27r/ da.Wp{a)Wp{ot). By linearity we can extend this 
to any mixed state p = ll/'XV'l- 


3. Symplectic covariance 


Qj = G^*Wp. (AlO) 


Under a linear symplectic transformation Y“j,, an ar¬ 
bitrary symplectic tensor with both covariant 

(lower) indices and contravariant (upper) indices trans¬ 
forms as 


Likewise, when the generalized Glauber P function exists 
it satisfies 

Wp = G^*pJ. (All) 


/^ab-" 

^ gh- 




.ya Y 


h' 


(A5) 


An equation with all upper and lower indices contracted 
together is automatically invariant under a (linear) sym¬ 
plectic transformation. For example, 

AcB‘^ AaY“ (r“i)=fcB^ = (A6) 


Appendix B: Damped harmonic oscillator example 

Here we calculate the characteristic time T and the 
covariance matrix V describing the POVM-and-prepare 
channel associated with the common case of a lightly 
damped harmonic oscillator. For a concrete example, 
consider the center-of-mass motion of the single cooled 
ion in a Paul trap described by Bushev et al. [SH]. In the 
absence of feedback, the master equation is 

dtp = -tiy[a^a, p] + r(iV + l)Va[p] + P/VDat [p] (Bl) 


Dynamical equations ([^ and (A4) exhibit such manifest 
symplectic covariance. 


The operations of lowering and raising indices with 
are compatible with an overall linear symplectic transfor¬ 
mations because they are characterized by their preserva¬ 
tion of the symplectic form, = Y® a- (The Wigner 
function is a scalar phase-space density which transform 
trivially under a linear symplectic transformation.) 


where a = (x^/m^'+ ipjsJmv)l\/2 is the harmonic trap 
lowering operator, m is the ion mass, v = 1 MHz is the 
trap frequency, P = 400 Hz is the laser cooling rate, N = 
{d)d) sa 17 is the steady-state occupation number, and 
the superoperator is defined by 

T>c[p] = cpc'l'- i{c'''c,p}. (B2) 
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Choosing length units such that -Jmvjh = 1, this can be 
brought into our standard form Q with 

11 = vl, D = r(Af+1/2)1, 7 = r/2. (B3) 


(Note that |D| > 7 ^ since N > 0.) Using Eq. ([^ and the 
defining equation |Ct| = (l+e“^'’'^)^/4 for T, we compute 


The dynamics (Bl) thus generate the evolution 


Po Pt = (d>j o Ur) [po] 


with 


^ _ ^-arccoth(2A^+l) 

_ _ ^Te _ ( cos vT sin z/T \ 

^ \ - sin uT cos i^T j ' 


(B7) 


(B 8 ) 

(B9) 


Ct = (7V+l/2)(l-e-‘r)I, 
2 

T = — arccoth {2N + 1), 

V = I. 


(B4) This is pure harmonic evolution for a time T followed 
by a POVM measurement in the Gaussian basis of co- 
(B5) herent states |l;a), where the corresponding prepared 
states are contracted toward the origin with proportion- 
ality constant s < 1 . 




